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Abstract 

We study the Cauchy problem for the Korteweg de Vries (KdV) equation with small dispersion 
and with monotonically increasing initial data using the Riemann-Hilbert (RH) approach. The 
solution of the Cauchy problem, in the zero dispersion limit, is obtained using the steepest de- 
scent method for oscillatory Riemann-Hilbert problems. The asymptotic solution is completely 
described by a scalar function Q that satisfies a scalar RH problem and a set of algebraic equations 
constrained by algebraic inequalities. The scalar function Q is equivalent to the solution of the 
Lax-Levermore maximization problem. The solution of the set of algebraic equations satisfies the 
Whitham equations. We show that the scalar function Q and the Lax-Levermore maximizer can be 
expressed as the solution of a linear overdetermined system of equations of Euler-Poisson-Darboux 
type. We also show that the set of algebraic equations and algebraic inequalities can be expressed 
in terms of solutions of a different set of linear overdetermined systems of equations of Euler- 
Poisson-Darboux type. Furthermore we show that the set of algebraic equations is equivalent to 
the classical solution of the Whitham equations expressed by the hodograph transformation. 

1 Introduction 

The Cauchy problem for the Korteweg de Vries (KdV) equation 

u t - 6uu x + e 2 u xxx = 0, u(x, 0) = u (x), e>0 (1.1) 
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in the zero-dispersion limit has been widely studied. The physical interest in this limit is due to the fact 
that it describes the phenomenon of shock waves in dissipationsless dispersive media. Dispersive shock 
waves are characterized by the appearance of rapid modulated oscillations. Gurevich and Pitevskii jj] 
suggested that these oscillations could be modeled by the solution of the one-phase Whitham equations 
Q . The multiphase or g-phase Whitham equations were derived by Flaschka, Forest and Mc Laughlin 
H . Lax and Levermore Q] rigorously showed that the multiphase Whitham equations appear in the 
zero dispersion limit of the Cauchy problem for the KdV equation with asymptotically reflectionless 
initial data. Later Venakides B considered a wider class of initial data. Lax and Levermore developed 
their theory in the frame of the zero-dispersion asymptotics for the solution of the inverse scattering 
problem of KdV. They showed that the principal term of the relevant asymptotics is given by the g- 
phase solution Js) of the KdV equation with the wave parameters depending on the functions u% (x,t) > 
■ • • > U2 g +i(x, t) which satisfy the g-phase Whitham equations: 

du du ' 

— -Vi(ui,U2,...,U2g + i)-^ =0, x,t,u l eJR, i = 1, ...,2g+ 1, g > 0. (1.2) 

For g > the speeds Vi{u\, U2, ■ ■ ■ , U2 S +i), i — 1, 2, . . . , 2g + 1, depend through ux,...,U2 g +i on 
complete hyperelliptic integrals of genus g. For j = 0we define u\ = u and the zero-phase Whitham 
equation reads 

du „ du .„ „, 

m- 6u d- x = ° (L3) 



and can be integrated by the method of characteristics. The formal integrability of equations ( |l.2[ ) 
for g > was obtained by Tsarev Q using the geometric-Hamiltonian structure || of the Whitham 
equations. Namely he proved that if the functions u>i — Wi(ui, U2, ■ ■ ■ , U2 S +i), i = 1, . . . ,2g + 1, solve 
the linear over-determined system of equations 

dwi 1 dvi 



duj Vi — Vj duj 1 



], i,j = 1,2,..., 2.9 + 1, i^j, (1.4) 



where V{ = Vi(u\, U2, ■ ■ ■ , U2 g +i), i = 1, • ■ .,2g + 1, are the speeds in (1.2), then the solution u(x,i) 



(ui(x,t),v,2(x,t), . . . ,U2 g +i{x,t)) of the so called hodograph transformation 

x = —Vi(u) t + Wi{u) i = 1, . . . , 2g + 1 , (1-5) 



satisfies equations (1.2). Conversely, any solution (ui(x, t),U2(x, t), . . . , v,2 g +i(x, i)) of (|l.2| ) can be 
obtained in this way in the neighborhood of (xo,to) at which it^'s are not vanishing. The general 
solution of the Tsarev equations was obtained in |j| , |hJ , [0 for monotonically increasing initial data. 
The key step introduced in || . |llj was to reduce the solution of the Tsarev system to the solution 
of linear-overdetermined systems of Euler-Poisson Darboux type for some functions = (u) , k = 
l,...,g, namely 

d d , . d 2 

T^—qk - T^—qk = 2{u t - Uj) gfc, i t i, l, J = l,...,2g+ 1. 

dui duj duidu-j 



2 



In this paper we show that these functions play a role in the Lax-Levermore maximization problem. 
We study the small dispersion limit of KdV for monotonically increasing analytic initial data bounded 
at infinity. 

We use a different approach to the small dispersion limit obtained in by Deift, Venakides and 
Zhou. These authors first used the formulation of the Cauchy problem for KdV as a Riemann-Hilbert 
(RH) problem |Q (see also |l4|]). Then, for computing the small e-asymptotics, they used the steepest 
descent method for oscillatory Riemann-Hilbert problems introduced in [ fl6) . Their procedure leads 
to a scalar RH problem for a certain phase function Q that turns out to be equivalent to the solution 
of the leading order variational problem in the Lax-Levermore theory. Furthermore they reduce the 
initial value problem for the Whitham equations to solving a set of algebraic equations constrained 
by algebraic inequalities. Existence and uniqueness of the initial value problem for the Whitham 
equations follows from the existence and uniqueness of the solution of the variational problem in the 
Lax-Levermore theory. 

The RH problem for the scalar function Q is well defined even for smooth initial data. 
We show that, for smooth monotonically increasing initial data bounded at infinity, the function Q 
and the Lax-Levermore maximizer can be expressed as the solution of a linear overdetermined system 
of Euler-Poisson-Darboux type. In the same way, the set of algebraic equations obtained through the 
Deift, Venakides and Zhou approach can be expressed as the solution of a set of linear overdetermined 
systems of equations of Euler-Poisson Darboux type. We also show that this set of algebraic equations 
is equivalent to the set of algebraic equations defined by the hodograph transformation (1.5). 
The advantage of this representation is clear when one tries to construct effectively the solution of 
the set of algebraic equations constrained by the algebraic inequalities. It is simpler to evaluate space 
derivatives and to estimate the sign of the quantities in the inequalities. Indeed this representation 
was used in |lO[],[{L7| to give an upper bound to the genus of the solution of the Whitham equations. 

This paper is organized as follows. 
Section 2 contains the definitions of the Abelian differentials on Riemann surface and the meromorphic 
analogue of the Cauchy kernel. 

In Section 3 we review the Riemann-Hilbert steepest descent method for the zero dispersion KdV 
equation with monotonically increasing initial data. 

We determine the scalar function Q associated to the RH problem and all its properties in section 
4. We show that the set of algebraic equations that we call moment conditions and normalization 
conditions are equivalent to the set of algebraic equations defined by the hodograph transformation. 
In Section 5 we show that the hodograph transformation obtained in Sec 4 is equivalent to the classical 
one provided in ]9| and We then show that the moment conditions, the normalization conditions 
and the Lax-Levermore maximizer can be expressed in terms of solutions of linear overdetermined 
systems of Euler-Poisson-Darboux type. We also derive in a simple way the equations which determine 
the phase transitions. 

Finally in section 6 we summarize the main results and draw our conclusions. 
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2 Riemann surfaces and Abelian differentials: notations and 
definitions 



Let 



2g+l 



S g := |P= (A, j/), y 2 = JJ (A-Uj)| , ui > m 2 > •• • > u 2g+l , Ui e H, (2.1) 

be the hyperelliptic Riemann surface of genus g > 0. We shall use the standard representation of S g 
as a two-sheeted covering of CP 1 with cuts along the intervals 



[«2*,«2fc-l], fc = l,...,flf + l, U2 9 +2 



-00 . 



(2.2) 



We choose the basis {aj, f3j} 9 j =1 of the homology group Hx(S g ) so that a>j lies fully on the upper sheet 
and encircles clockwise the interval [v,2j, v,2j-i], j — 1, ■ ■ ■ , <?, while /3j emerges on the upper sheet on 
the cut [v,2j, v,2j—i], passes anti-clockwise to the lower sheet trough the cut (— oo,U2 g +i] and returns 
to the initial point through the lower sheet. 

The one-forms that are analytic on the closed Riemann surface S g except for a finite number of points 

are called Abelian differentials. 

We define on S g the following differentials (TsJ : 

1) The canonical basis of holomorphic one- forms or Abelian differentials of the first kind <j>\, 4>2 ■ • ■ <fig : 



+ X g-2 1 k + .,. +1 k 



d\ , k = 1 , 



■!)■ 



(2.3) 



The constants 7* are uniquely determined by the normalization conditions 

fo = Sjk, 3,k = l,...,9- (2.4) 

We remark that an holomorphic differential having all its a-periods equal to zero is identically zero 
@. 

2) The set af., k > 0, g > 0, of Abelian differentials of the second kind with a pole of order 2k + 2 at 
infinity, with asymptotic behavior 



and normalized by the condition 



We use the notation 



fc-i 



0{\~ 



dX for large A 



of = 0, j = l,...,g. 



(2.5) 



(2.6) 



<rg(A) = dp9(\) , 12oi(A) = dq 9 (X) g>0. 



(2.7) 
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In literature the differentials dp 9 (X) and dq 9 (X) are called quasi- momentum and quasi-energy respec- 



tively p). The explicit formula for the differentials a 9 k , k > 0, is given by the expression 
°iM = , P 9 (X) = A«+ fe + a^A^ 1 + a k 2 X 9+h - 2 ■ ■ ■ + a k g+k , 



where the coefficients af = af(u) , u = (u%, U2, ■ ■ ■ , u^g+i), % = 1, ...,<? + k, are uniquely determined 



by and (2.6) 



3) The Abelian differential of the third kind oj qqo (X) with first order poles at the points Q = (q, y(q)) 
and Qo = (qo,y(qoj) with residues ±1 respectively. Its periods are normalized by the relation 

uj qqo (X) = 0, j = l,...,g. (2.9) 

In the following we mainly use the normalized differential wf (A) which has simple poles at the points 
Q ± (z) — (z,±y(zj) with residue ±1 respectively. 

The differential ui 9 z (A) is explicitly given by the expression 

dX y(z) A /" dt y(z) 
wf(A) = -TTT-- Vft A -7TTT , (2-10) 

where </f>fc(A), = 1, . . . , g, is the normalized basis of holomorphic differentials. By construction 

wf(A) = 0, j = l,..., ff . (2.11) 

The differential w|(A) can also be written in the form 



where 

dry 



iV^,u) = ^)^ 7 i / 
.7=1 ^ 



y(»?)(»7 - 2) 



(2.13) 



wf (A) as a function of z, is an Abelian integral. The periods of this integral are obtained from the 
relations fll9l 



4K(A)] = 0, / d,[wf(A)]=47ri0 fc (A), ./ 1- •••'/• (2.14) 

The differential uj 9 (A) satisfies the property jl9| 

d,wf (A) = (2.15) 
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In the following we will use the single value restriction of u> 9 z (A) determined by the conditions 

wf(A)U=„ 29+1 =0 



and by choosing y{z) to be analytic off the cuts (2.2) and real positive z > u\. We will still denote this 
single value restriction with (A). We remark that wf (A) is a meromorphic analogue of the Cauchy 



kernel on the Riemann surface S g |lj 



The next proposition is also important for our subsequent considerations. 



Proposition 2.1 The Abelian differentials of the second kind cr^(A), k > 0, defined in ([£ 
satisfy the relations 



4(A) 



- Res 

2 z— oo 



wf(A)z fe ~^z 



f 



-d\ Res 



(2.16) 



2fc + 1 Z = 00 

where (A) /ias Seen defined in uJ^(z) is the normalized Abelian differential of the third kind 

with simple poles at the points Q ± (A) = (A, ±j/(A)) wii/i residue ±1 respectively and d\ denotes 
differentiation with respect to A. 



3 Riemann-Hilbert steepest descent method for the zero dis- 
persion KdV equation with monotonically increasing initial 
data 

Following jl8| we reformulate the inverse scattering for the KdV equation as a RH problem. We 
consider monotonically increasing analytic initial data Uq(x) bounded at infinity. For convenience wc 
assume 



lim Uq(x) = 0, lim Uq(x) = 1. 



x — > — oo 



x — > + oo 



We suppose 



u (x){l + \x\ 2+s )dx < oo, 



(1 - u (x))(l + \x\ 2+s )dx < oo 



(3.1) 



for all finite c and 5 > 0. 

Let r(A;e), A > 0, be the reflection coefficient from the left of the Schodinger equation —e 2 f x 
u (x)f = A/. Define the matrix fTl, |l| 



-f e - 2la /^ 



v(X,e) = < 



re 



2ia/e 



1 



I - |r| 2 -f e - 2 ' iQ / £ 



A < 

< A < 1, 

A > 1, 
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where a\ = ^ and a = 4t\i + xX? . The goal is to find a row vector valued function m(A) = 

m(X;x,t, e) = (mi, 77x2) analytic for complex A off the real axis, satisfying the jump and asymptotic 
conditions 

m + (A; x, i, e) =rn_(A; i, e)f(A; a;, i, e) 
m(A; 2;, f, e) — >(1, 1) as A — > 00, 

where m± = lim,5^o m (A ± i8;x, t, e). The RH problem of finding the matrix m(A) given f(A) has a 
unique solution in the space (1,1)+ L 2 (dX?). The solution of the Cauchy problem (1.1) is given by 

u(x, t, e) = —2ied x mn(x, t, e) (3-2) 

where 

mx(X;x,t,e) = 1 + + 0(A~ 1 ), A -> 00, (3.3) 



see again (|18|, (14[). In this paper we study the Cauchy problem (1.1) in the limit e — > 0, the so called 
zero-dispersion limit of the KdV equation. We use the WKB approximation with one turning point 
to calculate the reflection coefficient |2Cj 

r(A;e)c-ze- 2 ^ A )/ e X[o,i](A) 

r (A) 1 

p(\)~\zx(\)— / [A 2 — (A — uo(x))?]dx, 
J — 00 

where the quantity x(X) is defined by the relation u (x(X)) — X. As usual X[o,i]W denotes the 
characteristic function of the interval [0, 1]. From the above considerations the jump matrix reduces 
to the identity matrix for A > 1 and our RH problem is reduced to the interval (—00, 1]. The quantity 
p(X) can be expressed also in the form 

(3.4) 

2 Jo vx-y 

where f{u)\t=o is the inverse function of the initial data uq(x). In the following we identify r with its 
WKB approximation. 

Following the procedure in Jl2[ , we introduce a change of the dependent variable m 

M (A) = m(A)e e(A)CT3/£ (3.5) 
where the scalar function G(X) = G(X; x, t) is analytic in A off the line (—00, 1] and satisfies G+(X) + 



£_(A) = 0, for A G (-00, 0) and Q(X) -> as A -> 00. From (gj) and (|3j) we obtain 

u(x,t) = -2iea 8 Mn(aj,t) - 23 x £i (x, t), 
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where£(A) = Qi(x, t)/y/\ + 0(1/ X) and M\(\;x,t) = 1 + M 11 (x 1 t)/\i. The RH problem in the new 
variable becomes M+(A) = M_(A)Y(A), A S (0, 1), where 



T = 



-ie~ th/e 



and h = Q+ + Q- — 2p + 2a. For computing the small e asymptotics we follow the technique in |L2 |. 
The interval < A < 1 is partitioned into finitely many intervals Ij = (u2j,U2j—i), j = 1, . . . ,.<? + 1 
and = U2g+2 < U2 g +\ < ■ ■ ■ < U2 < t«i < 1, g > 0. Using the steepest descent method |0|,[[L6|, the 
jump matrix T can be reduced to one of the two forms with exponentially small errors as e \ 0, 



oO) f °, / 16 \ Aeu'^-, 



60) °), Ae^iAuf^ 1 /, 



(3.6) 



The function Q satisfies the following conditions: 

<y<3+ < 0, and h'(X) = 0, thus G+ + Q- - 2p + 2a = -fi,-, A G U^ 1 /,-, (3.7) 

where Qj is some constant of integration; 



0+ - Q- = and fc'(A) > 0, A e (0, 1)\ \jSl\ (u 2j , u 2j -x) 



(3.8) 



The remaining RH problem is the following 

M + = M_(j 1 , A e (-oo,0), 



(3.9) 



where 03 



'1 
,0 -1 



To complete the solution of the initial value problem (1.1) we have to determine the phase function 
Q, the intervals Ij, j = 1, . . . ,g + 1, the values flj. Once all these steps have been made, we have to 
solve the RH problem (3.9), so that the solution of the initial value problem (1.1) can be expressed 

by H, f§ 



2g+l 



i(x, t,e)=J2 «j + 2a? - 2e 2 ^ log0(Q/(27re)), 
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where a® has been defined in (2.8) and il = (fii, f^, ■ ■ • , The theta function is defined by 



9(z,B)= 



cl 2ni('m,z)-\-7ri(m,B'm) 



meZs 



where B is the period matrix of the holomorphic differentials (2.2), namely B t j = j„ i, j = 1, . . . , g. 



4 Determination of Q{\) 

We observe that for each fix x and t the function G'(X) satisfies the following RH problem 

g' + + g'_=o, a € (-00,0) 

Q' + -QL=0, X> landAe (0, 1) - £ ff , (4.10) 
G' + + Q'_ - 2p' + 2a' =0, XeC g , 

where C g — U^Jij. We call the intervals Ij bands, while the intervals (u2j+i, «2j), j — 0,...,<?, 
mo = 1, are called gaps. We add the requirement that 

(\/A^'(A))±, are continuous functions for real A. (4-11) 

It follows that Q'(\)± are continuous functions for A € (0, 1). 

We recall that the condition G(X) = 0(X~i) for large A, implies 

g'{\) = O(X-i). (4.12) 



We observe that the condition (3.8) implies the following normalization condition 

g'(X)d\ = 0, j = l,..., g, (4.13) 



where ctj is any clockwise close loop around the cut Ij. When the loops ctj collapse to the interval Ij 
described twice, the conditions ( 4.13| ) become 



/ (g'+-gL)dx = o, ] = 



(4.14) 



The solution of the RH problem ( 4. 10] ) that satisfies ( |4.1l| ) is given by the integral J22 



r , m ^(A) f (2p'(z) - 2a'{z))dz 



where 

2g+l 



2 



Y[(X- Uj ). (4.16) 

3=1 
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We choose y(X) to be analytic off the intervals (— oo, 0] U £ g , real and positive for A > u\ and we 
denote y + (X) the boundary value from above the cut (—00. 0] U C g . 

In order for ( 4.1ii ) to satisfy ( 4.12| ) we must impose the following moment conditions 



M1^a = o, k = 0,. 
y + W 



(4.17) 



Furthermore we must impose the normalization conditions (4.13). We observe that ( 4.13| ) and (4.17) 
represent a system of 2g + 1 algebraic equations which, in principle, determines the end-points 
Mi, iti, . . . , U2 g +i of the intervals Ij, j = 1, . . . , g + 1. 
For 5 = the moment condition reduces to the form 

p'(X)-a'(X)) 



2 
ni 



Vx^ 



-d\ = f(u) — 6tu — x = 0, 



which is the solution of the zero-phase equation (1.3). 
We observe that 



Q' + (z Q )-g'_{z ) = g{z Q ) = 



y(z ) f (2p'(z) - 2a'(z))dz 



Til 



(z - z )y+(z) 



and 



0+(zb)-a_(3b) = o, Ae (0,l)\uf+i (u 2j ,u 2j -i). 

Therefore because of the assumption of the continuity of Q'± on (0, 1), the following relation must be 
satisfied 



lim g'(z ) = 0, i = l,...,2g + l, 



(4.18) 



where zq € (ui + i,Ui) for i odd and zo € (^i-ij^i) for « even. 

We consider analytic initial data or smooth initial data in (0, 1) that satisfies ( |3.l|) , so that the function 
p'(X) is Holder continuous in subsets of (0, 1), namely 

\p'(Xi) - P '(X 2 )\ < c\Xx - X 2 \ s , VA 1; A 2 e J 

where the constant c > 0, < <5 < 1, and J is some open subset of (0, 1). 

Lemma 4.1 Jj||/ // the function p'(X) is Holder continuous near and at Ui, i = 1, . . . , 2g + 1, then 

lim g'{z ) = 0, i = l,...,2g + l, 

where zq £ (ui + \,Ui) for i odd and zq G (ui^i,Ui) for i even. 
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The above lemma guarantees the consistency of the assumptions (4.11) and (4.18). 
These considerations suggests to build a second solution Q' (A) of the RH problem (4.1C) requiring 
that G'(X) satisfies ( 4.12 ) and ( 1.14 ) while we impose the continuity ( 4.11 ) as constraint. The solution 
Q' (A) is given by the expression 



G'W = 



(2ni)y(X) 



2y+{z){p'(z)-2a'(z))dz Q fl _i(A) 
z-X y(X) ■ 



(4.19) 



The polynomial Q g -i(X) has degree g — 1 and its coefficients are uniquely determined from (4.14) or 
(4.13). It is easy to verify that 



Qg-iW 



1 9 

— y 



2ni ^-^ dX 

k=l 



dz(2p'(z) - 2a'{z)) 



y + {z)drj 



(4.20) 



where 4>k{X), k = 1, ... g is the basis of holomorphic differential defined in fl2.3p . 

Lemma 4.2 The function Q'(X) satisfies the conditions (4-H) and 

The continuity on the function Q 1 (A) is obtained imposing that the end points u\, 
according to the equations 

S'(A) A=Ui = 0, « = 1,...,2<7 + 1. 



,u 2g +i evolve 



(4.21) 



The next theorem establishes the equivalence between the two different solutions Q' (A) and Q' (A) of 
the RH ( p~To| ). 



Theorem 4.3 The function function G'(X) and the set of algebraic equations (1-21) is equivalent to 
the function G'{X) and the set of algebraic equations (4-l r i) and {{.th). The equivalence is established 
for any C°° initial data satisfying (3.1). 

Proof: We write G'(X) in the form 



G'W = 



1 



Using the identity 



y 2 (z)(2p'(z) - 2a'(z))dz 
(2iri)y+(z)(z-X) 



E 

i=i 



dX 



dz 



{2p'{z) - 2a'(z)) 
{2m)y+(z) 



- X 



k-l 

^va*- 1 - 

3=0 



y 2 (z)dn 

y{v){z-v)' 

(4.22) 



(4.23) 



we obtain 



y\z) y\X) 



2tf 



A 



A 



j=0 



(4.24) 



k=0 
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where the s^'s are the symmetric function in the variables Ux, 112, ■ . ■ U2 g +i, namely sq — 1, s\ 
Y^k=i u ki s 2 — J2k<j u kUj and so on. Using ( 4.24 ) we can rewrite ( 4.22 ) in the form 



y{\) f (2p'(z)-2a'(z))dz 



2ixi 



y+(z)(z-X) 



E 

1=1 



&(A) f (2p'(z)-2a'(z)) 



{2ni)d\ 



dz 



z — rj 



dr] 



1 2g k , 



(2m)y(\) 



k=0 
2g 



3=0 
k 



(2p'(z) - 2a'(z))z j 



dz 



(4.25) 



1=1 fc=0 j=Q J ^s 



(2p'(z) - 2a'{z))zi 



dz 



k-j 



dr/. 



Imposing (4.21), we can see that the first term in (4.25) is automatically zero at the branch points by 

) imr 



lemma 4J. Therefore, using ( 4.25 ), the equations ( 4.21 ) imply 

n n— 1 

/ I ■'/>' I - ) — ''it' I - I I / ni 1: I 

-dr\ 



E E 7fA m / 



(2p'(z) - 2a'(z)) 



1=1 m=0 
2ff 



y+(z) 



dz 



Z — f] 



E(-) fc ^E^7 

f n a — n J C 



(2p'(z) - 2a'(z))z j 



k=0 

g 9-1 



29 



y + (z) 



dz 



EE^Er^W 

l—l m=a k=0 j=0 C 9 



{2p'{z) - 2a' '{z))z 3 

yW) 



dz 



r) k - j 

yiv) 



dr) 



0, 



for i — 1, . . . , 2g + 1. The above quantity is a polynomial in the A variable of degree 2g that must 
have 2g + 1 zeros, therefore it is identically zero. From the coefficients of degree g to 2g we get the 
moments conditions (4.17). For the coefficients of degree m, m = 0, . . . , g — 1 we get the relations 



n 

1=1 



(2p'(z)-2a'(z)) 

y + (z) 



dz 



yiji) 

z — rj 



dr\ 



2fl 



EtT E (-)W-* E 



i=l fc=g+l 



" ! ~ k=m+g+l 

(2p'(z) - 2a'(z))z j 



29 r 

E / 



(2^(2) - 2a'(z ))z fc 

y+(*) 



k—ra 



-dz 



y + {z) 



dz 



y(v) 



drj = 0, to = 0, . . . , g — 1. 



Using (2.4) the above relation simplifies to the form 



(2p'(z) - 2a'(z)) 

y + {z) 



dz 



z — rj 



drj = 0, to = 0, . . . , g — 1. 



(4.26) 



Because the matrix {7™}, I = 1, . . . , g, m = 0, . . . , g — 1 is invertible, the relation ( 4.26| ) is equivalent 
to 



L 



(2p'(z) - 2a'(z)) 

y + (z) 



dz 



y(v) 

h*-V 



dr] = 0, 1 = 1, 



,9 



(4.27) 
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which coincides with the normalization conditions (4.13). Therefore on the solution of (4.17) and 
(4.13) the function Q 1 (A) reads 



G'W = 



y(X) f (2p'(z) - 2a'(z))dz 



2ni 



y+(z)(z-X) 



= Q'{x). 



(4.28) 



□ 

In the same way we can prove that £7'(A), the moment conditions (4.17) and the normalization condi- 
tions (4.13) are equivalent to G'(X) and the equations (4.21). We remark that the equivalence between 
the two solutions does not depend on the fact that the problem has a priori a unique solution as 
follows from the Lax-Levermore theory but it derives only on the structure of the RH problem for the 
scalar function Q . 

We observe that using oj 9 z (A), the meromorphic analogue of the Cauchy kernel defined in ( 2.10 ), we 
can write Q 1 (X)dX in the form 



g'(x)dx = ~xdf> 9 {x) - tdq 9 {x) + n»(A), 



(4.29) 



where 



dp 9 (X) =- 



2iri 



<j*(X)z- 



and 



n»(A) = ~ [ 2"l(X)p'(z). 

27rl J Co 



Here and below the integrals in the z variable are taken on the upper side of C g . 
Using the above representation we compute the constants defined in (3.7) 



g'(x)dx 



0k 



2-ki 



Integrating by parts the above identity and using fl2J5| ) and ( |2li| ) we obtain 



tt k = 4 



(j) k {z){{p{z) - a(z))dz, k = l,.. 



From the above it follows that fi s +i = 0. 
The following theorem due to Krichever 



(4.30) 



(4.31) 



(4.32) 



connects the solution of the set of algebraic equations 



(4.17) and (4.13) or (4.21), to a solution of the Whitham equations. 

Theorem 4.4 J§q/ Let us suppose that the Ui 's depend on x and t in such a way that the conditions 



(4- 21) are fulfilled. Then m = Ui(x,t), i = 1, . . . , 2g + 1, satisfies the Whitham equations 

= 1,.. .2.9 + 1, 



d d 



(4.33) 
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where 



v t (u) 



dq3(\) 
dp9(\) 



■ 2.9 + 1 



(4.34) 



A— Ui 



and the differentials dp 9 (X) and dq 9 (X) have been defined in ({.31 ). 



We remark that from proposition (2.1) the following identity is easily verified 



Vi{u) 



dq 9 (X) 



dpa(X) 



dq 9 {X) 



X—Ui 



d P 9(X) 



i = l,...,2 ff + l, 



(4.35) 



A— Ui 



where dp 9 (X) and dq 9 (X) have been defined in (2.7). The second expression of the Vi(u)'s in (4.35) is 
the classical formula for the speeds of the Whitham equations obtained in m . 



We can write the algebraic equations ( |4.2lD in the form of the so called hodograph transformation 
introduced by Tsarev |7|: 



where 



Wi{u) 



■Vi(u)t + Wi{u), i = 1, . . . ,2g + 1, 
n 9 {X) 



dp9{X) 



1,...,2«? + 1 



(4.36) 



(4.37) 



A— Ui 



and fl 9 (X) has been defined in (4.31). As a consequence of theorem 4.3, the set of algebraic equations 



( 4.17 ) and ( 4.13 ) is equivalent to the hodograph transformation ( 4.36 ) 



In the next section we will show that the u>i(u)'s defined in (4.37) coincide with the classical formulas 
provided in ||] or |l(J ■ 



5 Solution of the Tsarev system and linear overdetermined 
system of Euler-Poisson-Darboux type 

We first define the Cauchy problem for the Whitham equations. The initial value problem consists of 
the following. We consider the evolution on the x — u plane of the initial curve u(x,t = 0) = uq(x) 



according to the zero-phase equation (1.3). The solution u(x, t) of (1.3), with the initial data uq(x), 
is given by the characteristic equation 



x = —6tu + f(u) 



(5.1) 



where f(u)\t=o is the inverse function of the initial data uq(x). The solution u(x,f) in (5.1) is globally 
well defined only for < t < t , where to = \ min„ e ]R,[/ / (u)] is the time of gradient catastrophe of 



( p.lD . Near the point of gradient catastrophe and for a short time t > t , the evolving curve is given 
by a multivalued function with three branches u±(x,t) > 1*2(2, i) > u^{x,t), which evolve according 
to the one-phase Whitham equations. 



14 



Outside the multivalued region the solution is given by the zero-phase solution u(x, t) defined 
in (5.1). On the phase transition boundary the zero-phase solution and the one-phase solution are 
attached C 1 -smoothly. 

Since the Whitham equations are hyperbolic pj| , other points of gradient catastrophe can appear in 
the branches 1*1(2;, > U2(x,t) > u 3 (x,t) themselves or in u(x,t). 

In general, for t > to, the evolving curve is given by a multivalued function with and odd number 
of branches u\{x,t) > 1*2(2;,/;) > ■ ■ ■ > U2 g +\{x,t), g > 0. These branches evolve according to the 
g-phasc Whitham equations. The <?-phase solutions for different g must be glued together in order to 
produce a C 1 -smooth curve in the (x, it) plane evolving smoothly with t. The initial value problem 
of the Whitham equations is to determine, for almost alH > and x, the phase g(x,t) > and the 
corresponding branches U\(x, t) > 1*2(2:, £)>•••> 1*23+1(2;, t) from the initial data x = f(u)\ t=0 . 
The solution of the Whitham equations for a given g is obtained by the so called hodograph transfor- 
mation introduced by Tsarev 0. 

Theorem 5.1 Ifwi(u), u= (ui, u 2 , ■ ■ . , 1*23+1); solves the linear over- determined system 
dwi 1 dvi 

— = — [wi-Wjl, i,j = 1,2,..., 2# + l, i^j, (5.2) 

OUj Vi — Vj OUj 

then the solution (ui(x, t), 1(2(2:, t), . . . , 1*23+1(2;, t)) of the hodograph transformation 

x = —Vi(il) t + Wi(u) i = 1, . . . , 2g + 1 , (5-3) 

satisfies system ( fi.j}j . Conversely, any solution (u\, 1*2, . . . , 1*23+1) of can be obtained in this way 
in a neighborhood (xo,to) where the Ui x 's are not vanishing.. 

To guarantee that the g— phase solutions for different g are attached continuously, the following natural 
boundary conditions must be imposed on Wi(ui, 1x2, . . . , 1*23+1), i = 1, . . • , 2g + 1, g > 0. 
When it; = it;+i, 1 < I < 2g, 

wf(ui,. ..,Ul,Ul,... ,U2 ff +l) = V)f +1 (ui, ...,Ul,Ul,... ,U2 ff +l) ( 5 -4) 

and for 1 < i < 2g + 1, i^ 1,1 + 1 

wf(ui, . . . ,m,ui, . . . ,U2 9 +i) = wf~ l {ui, . . . ,ui,ui, . . . ,it 2g +i). (5.5) 

The superscript g and g — 1 in the w^s specify the corresponding genus and the hat denotes the 
variable that have been dropped. When g = 1 and 1*2 = 113 we have that 

ifi (1*1,1*3,1*3) = /(tti) 

W2(Ul,U3,U 3 ) = W 3 (l*i, 1*3, 1*3) , 

where f(u) is the initial data. Similar conditions hold true when 1*1 =1*2, namely 

11)3(1*1,1*1,1*3) = /(i* 3 ) 

Wi{Ui,Ui,U 3 ) = W2(Ui,Ui,U 3 ). 
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We remark that the Ui(tt)'s satisfy the boundary conditions (5.4-5.5) and for g = 1 we have 

vi(ui,u 3 ,u 3 ) = -6ui, u 3 (ui,ui,u 3 ) = -6u 3 . 

The solution of the boundary value problem ( |5.2| ) , ( [5~4| - |5.7| ) has been obtained in for any smooth 
monotonically increasing initial data. 



Theorem 5.2 Let be f(u) the inverse function of the smooth initial data uq(x,0). If the function 
Qk = Qk(ui, v-2, ■ ■ ■ ! U2g+i)> 1 < k < g, is the symmetric solution of the linear over-determined system 

d 2 q k (u) dq k (u) dq k (u) 



2(ui - Uj) 



duiduj 



On, 



i^j, i,j = l,...,2g+l, g>0 



q k (u,u, ...,u) = F k (u) 

2g+l 



(5.8) 



F k (u) 



2(9- 



1 d g ~ 



( u g-hf{*-V{u)\ , 



(2g-l)ir du3' k 

with the ordering 1 > u\ > u 2 > ■ ■ ■ > it 2 g+i > 0, then Wi(u), i = 1, . . . , 2g + 1, defined by 



Wi(u) 



2g+l g k 

2d Ui q g (u) Yl i u i ~ u «) + X! 5Z( 2n ~ 1 )r/ s -n-P^_i('"0 

n— l.n^i fc— 1 n—1 

solves the boundary value problem (5.£), (5.4-5//). Conversely every solution of (5A), (5.4-5.%) 
be obtained in this way. 



(5.9) 



can 



In ( |5.9| ) the polynomials P%'s have been defined in fl2.8| ) and the IV s are the coefficient of the expansion 
for £ — > 00 of 



(5.10) 



The existence and uniqueness of the solution of the boundary value problem ( |5.8| ) has been proved 
in IT1. 



Theorem 5.3 |_?J|/ TTie solution of the boundary value problem ( |5. j| j is unique, symmetric with respect 
to the variables U\,U2, ■ . ■ , U2 g +i and reads 

-1 /•! /•! 



<7fc(0 =7^ , , 



dzidz 2 ...dz 2 9(l + Z29) 9 ^l + ^g-l) 9 * • • • (1 + Z 3 ) * (1+ Zl) 



F fe (±^(. . . (l±a(i±£i Ul + l^iua) + i^ U3 ) 



^r^ 2s ) 



v /(l-z 1 )(l-z 2 )...(l-z 2ff ) 



(5.11) 
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where C — Yli=i Of and 



C, 



1 + 

1 \/l - i" 



(5.12) 



The functions F k (u) and the solutions q k (u), k = 1, . . . , g, of the boundary value problem (5.8) satisfy 
the following relations: 

2.9 + 1 



duF k (u) 



-F k+ i(u) + ud u F k+1 {u), k = l,...,g-l, g > 0, 



(5.13) 



d Ui qk{u) = -q k +i{u) + Uid Ui q k +i(u) i = l ,25 + 1, k=l,...,g-l g > 0. 



Lemma 5.4 JldjJ The solution of the Whitham equations described by (5.5) where the Wi(u)'s 
given by (pM) is C 1 -smooth on the phase transition boundaries. 



The next theorem shows the equivalence between the hodograph transformation defined in (4.36) 
and the one define in (5.3) 



Theorem 5.5 For any smooth monotonically increasing initial data satisfying ( 3. 1 ), the following 
identity is satisfied 



w^u) = Wi(u), i = l,...2g+l, 



(5.14) 



where the Wi(u) 's are defined in ( 4-3V ) and the Wi(u) 's are defined in ( 5.L ). 



Therefore combining theorem 4.3 and and theorem |5.5| , we deduce that the hodograph transfor- 
mation (5.3) is equivalent to the set of algebraic equations (4.17) and (4.13). 

The next theorem shows that the hodograph transformation (5.3) can be written in a nice algebraic 
form. This is the first step to transform the moment conditions (4.17) and the normalization conditions 
(4.13) into a combination of solutions of linear-overdetermined systems of Euler-Poisson-Darboux type. 



Theorem 5.6 For g > the hodograph transformation \5.S\ ) where the Wi(u)'s are defined in ([^ 
is equivalent to the following set of2g + l algebraic equations 

23+1 

^ d U] q g - k {u) - kq g - k+1 (u) =0, k=0,...,g-2 

23+1 

d uM U ) ~ (-9 - - 6t = 0, (5.15) 

3=1 

23+1 2 9 +l 

2 2J Ujd Uj qi(u) + gi (u) - x - 6 1 ■ = 0, 

3=1 3=1 
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y(A)$(A;w)dA = 0, k=l,...,g 
where the function $(A; u) is given by the relation 



(5.16) 



2g+l 



(5.17) 



The function ^ 9 (X;u) satisfies the linear overdetermined system of Euler-Poisson-Darboux type 



^*.(A;iO-^(A;tO=2(« i -« J )^^(A;tO > i* i, i,i = l,..2 fl + l 
^-*»(A; fi) - 2^-**(A;iZ) = 2(A - Uj .)^A_*ff(A;a), j = 1, . . . 2 5 + 1 



SA 

*P(A;A,...,A) 



(5.18) 



2-" 



2g+l 



(29 + 1)!! 



/00(A) 



where /^0(A) is </ie gr£/i derivative of the smooth monotonically increasing initial data f{u). 
The above boundary value problem can be integrate in the form 

^ 9 (\-u)=-^J y -..y rfz 2 dz 2 ...dz 29+2 (l + z 2g+2 ) ff (l + z 2g+ i)^5...(i + Z3 )lx 

/to) ( l±^±a ( . . , (l±a(l^A + ^ Ul ) + l^ug) + -..) + ^fg^ugg+i) 
V(l- Z2 )(l-^)...(l-^ s+ 2) 



(5.19) 



where if = Yl ^j^i 0/ an d the Cj 's have been defined in (|5.12|) . The solution obtained is symmetric 
with respect to the variables u\, . . . , U2 g +\- A similar formula can be obtained for <E> 9 (A; u). 

In the next section we identify the system (5.15) with the moment conditions (4-17) and system 
( |5.16 ) with the normalization condition (4.13). 

Proof of Theorem |5.5| . 

We show that the Wi(u)'s defined in ( 4.37 ) satisfy the Tsarev system (5.2) and the boundary conditions 
(5.4-5.7). Therfore by theorem 5.2, we obtain Wi{u) = Wi{u). We follow the steps in J9[ , p3[] , 

The quantities Wi(ui, u 2 , . . . , U2 g +x), i = l,...2g + l, are well defined. Indeed let us write dp 9 (X) 

P 9 (v v) y 9 (\ il) 

and fi»(A) in the form dp 9 {\) = Vrf^ dA and fi9 ( A ) = ,A dX where 



y(A) 



X 9 {\U) = ; 



A — z 



(z)-^A^ [ N?(z)p>(z 



3 = 1 



(5.20) 



18 



and the N 9, s have been defined in ( |2.13 ). Here and below all the integrals are taken on the upper side 



of C g . Then 



Wi{Ul,U 2 , ■ ■ ■ ,u 2g +i) 



P9(u h u) 
1 



■KlP9{Ui, U) 



Jc g Ui-z j=1 J Cg 



(5.21) 



Next we show that the w^s satisfy the Tsarev system ( p.2| ). The differentials d Uj Q 9 (X) and d Uj dp 9 (X) 
are normalized Abelian differentials of the second kind with a pole at A = Uj of second order. Because 



d 



d 



of (4.37) the differential — — Q. 9 — Wj — — dp 9 (A) is holomorphic and 
N 1 ouj duj 



o =a, 

: 



d 



(n 9 (x)-wjdp 9 (x)) 
d 



n 9 



'd U j 



dp 9 (X), fc=l,. 



that follows from (2.11). Therefore 

(d U] n 9 (\)-w 3 d U] dp 9 (x)) = o 

because it is a holomorphic differential having all the a-periods equal to zero. 



From ( 5.22 ) we obtain 







d -, 



^-X 9 W-w j7 -P 9 (X) 



1 X g W - w,P 9 (X) 



duj A v ' 3 dun v ' 2 X — m 



, i= 1,2,..., 2.9+1. 



We use the above identity to evaluate djWi(u). From (5.21) and ( 5.23Q we obtain 
d _ ^ d x 9 (ui,u) . . 

— WAU) —— = , I + 7 

d Uj X 9 (ui ,u) -wj d Uj P 9 (im , u) „ <9„, P ff (itj , u) 

\-(Wj-Wi)- - 



P 9 (ui,u) 



P 9 (u l: u) 



which shows that 



_ d u P 9 {ui,u) l x 9 {u t ,u) - w 3 P 9 (ui,u) 

-(W-j — Wi) — = ^ 

P 9 {u z ,u) 2 (u i -u j )Ps(u i ,u) 

, _ _ . d Uj P 9 (ui ,u) Iwi- Wj 

--{Wj — Wi) — = -, 

V 3 ' P 9 (u h u} 2 Ui-uj 

1 dun d Uj P 9 (ui,U) 1 1 



Wi - wj duj P g {ui, u) 2ui — Uj 



(5.22) 



(5.23) 



(5.24) 



(5.25) 
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In particular the above argument also applied to dq 9 (X), dp 9 {\) and Vi(u) therefore we also have 



dvi d U] Pa{u u u) 



v l - vj duj P g (ui, u) 2 Ui - Uj ' 



which when combined with (5.25) proves the Tsarev relation for the w^s. 

Next we show that the functions Wi (u\ , u%, . . . , U2 S +i), i = 1, . . . , 2g+l, satisfy the boundary conditions 
([5.4H5.7D. 

In the following we use the superscript g to denote the corresponding genus of the quantities we 
are referring to. We need to consider the behavior of the Abelian differential wf (A) when two branch 
points become coincident. For the purpose let be ui = v + y/5, = v — y/5 where < 5 -C 1. 
The differential cj| (A) = lo 9 (A; u\, . . . , Ui—i, v — y/5, v + y/5, . . . , U2 g +i) has the following expansion for 
5^0 when I is odd (g| 

w°(\,8) = ^(A) + 6 -^-\v)d v ^ v {X) + 0(8 2 ), (5.26) 

where Wjj -1 (A) is the normalized Abelian differential of the third kind having first order poles at the 
points Q ± (z) = (z,±y(z)) with residue ±1 respectively and it is defined on the Riemann surface 

2g+i 

Sg-i-f = J] ( A -^)- ( 5 - 27 ) 

k— 1 



The differential 0(5 )/5 has poles at A = v of order at most 4 and zero residue. In formula (5.26) 
the quantity wf _1 («) = z v - U=t>- 

When 2 is even (A) has the following expansion |2j| for <5 — » 



^(A) - lo 9 z -\\) - -L^A) / Q '"V'W- (5-28) 



log 5 



The above expansion contains also terms of order y/~5/\og/5. Using ( 5.26| ) and (^.28| ) we can get the 
expansion of the differentials f2 9 (A) and dp 9 (X) when ui = v + y/5 and ui+\ = v — y/5, 1 < I < 2g. Let 
be C v the contour from v — y/5 to v + y/5 on the upper sheet of S g . When I is odd C g = C g -\ U C v 
where 

£g-l = [0, U2 3 +l] U • ■ ■ U [u l+3 , U l+2 ] U U Z _ 2 ] U ■ • • U [u 2 , 111] 



is the corresponding contour defined on the Riemann surface S g -\ of genus g — 1. From (5.26) we 
obtain the expansion of f2 9 (A) = fi 9 (A; ui, . . . , ui—\, v — y/5, v + y/5, . . . , U2g+i) as <5 — > 0, namely 

n 9 (x)^n 9 - 1 (x) + ^-d v uj 9 - 1 (x)n 9 - 1 (v)- - [ u g z (X) P '(z)dz + o(5 2 ), (5.29) 

2 m Jc v 
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where 



When / is even C g = £ g -i\C v where 



wl-\v) P '{z)dz. 



Cg-i = [0,«2 ff +i] U • • ■ U [ui +4 ,u l+3 ] U [u /+2 ,u/-i] U [u z _ 2 ,u z _ 3 ] U • • ■ U [u 2 ,Ui 



and from ( |5.28 ) we obtain the expansion of f2 9 (A) 



n 9 (\) 



uj 9 z (\)p'(z)dz + — 



a;f(A,J)p'(z)dz 



^(A) + — <-(A) 



-(«) i , 

n«- 1 («) + — / wf(A)p'(z)^. 



(5.30) 



The same expansions applies to dp 9 (X), therefore combining (5.29) and ( 5.30| ) we obtain 



fis(A) 



dp 9 (X) 



Of- 1 (A) 



r A=«i i dps- 1 (A) 

[ui=ui +1 =v\ 



, i ^ J.J + l, i = l,...2g + 1. 



(5.31) 



From the above identity it is clear that the boundary conditions ( |5.5| ) are satisfied. In order to evaluate 
(5.31) at the points A = v ± yS, we need to do some extra work. Let us defined the quantity 



X 9 (A)=y(A) 



dA ' 



(5.32) 



When = u + and = v — \/S, I odd, using (5.29) we obtain 



X 9 (A) =(A-«)x 9 - 1 (A) 



y(A) 



cf(A)p'(z)<fe 

9-1 



(5.33) 



+-tt»(d^(«) - (A - v) £ A"- 1 "^^- 1 ^)) + 0(5 2 ), 



fc=i 



where now 0(5 2 ) /5 2 is a polynomial in A. Using ( 5.26| ) we obtain the following expansion of the second 
term in the above equation 



y(A) 



u 9 {X)p'{z)dz =— 



1 f y + ^)j 



c v 



p (z)dz H : 

A — Z 7TZ 



(A-v) 



HI 



„ v/(z - u)2 - 5 



y r{z)p(z)dz 



fe=i 
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Therefore 



L0 9 (X)p'(z)dz 



1 

x=v±Vs m 

1 



~+( \ 1 1 ,Vz - v±V5 
V {z)p 0) . =dz 
G v \f z - v T VS 

y + {z)p'(z)dz 



7T 7 C „ ^/tf - (z - w) 2 



(5.34) 



m k=i c " 

= ±V*y + («y(«)+^ 1 (*), 
where lim^o /V1 ^ = 0. Combining ( 5.33| ) and (5.34) we obtain 

X 9 {v ±VS)= ±V5( X g - 1 (v) - y + (v)p'(v)) + TZ^S). 
In the same way we can get the expansion for P 9 (v ± V~5) , namely 



(5.35) 



P 9 {v ±VS) = ±VS(P 9 - 1 {v) 



2^ ' 



which, when combined with ( 5.35| ) gives 

x 9 (v± Vs) _ x a -Hv)-y + (v)p'(v) 



P9(v±VS) Pa- l {v)-y+{v)/{2^) 



0(6), 



(5.36) 



where lirn^o T^s(S) = 0. 

Using (5.3C) we obtain the following expansion of x 9 (A) when ui — v + and uj+i = v — VS, I even, 



X 9 (X)^(X-v)x 9 - 1 (X) + 



y(X) f W «(A) , 



TV l 

9-1 



p (z)dz 



+i^(y( u )-( A - u )E A9 " feiv r 1 («)) / Q W « 9 - 1 

lo S d fe= i •/«-(«) 

From the above and ( |5.34 ) we can evaluate 



(z). 



5-1 



X 9 (v ±VS)~ ±V5(x 9 - 1 (v) + y(v)p'(v)) - j^-Mv) ± V5J2 v'- 1 - k N«- 1 (v)) ^ ^ ^(z) 



k=l 



and 



P 9 (v ±V~S) ^ ±V~6(P 9 -\v) + ^y(v)) + ^(y(v) ±V~Sj2v 9 - 1 - k Nr 1 (v)) [° ^ , 
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Combining the above two expansions we obtain 



Pa 



(5.37) 



The relations ( 5.36 ) and ( 5.37 ) show that the boundary conditions (5.4) are satisfied. From ( 5.31 ) we 
obtain 



P ' {Z) Ay. 

iui(ui,?;,u) = -~- — - 



p°(ttl) f«l 

Jo 



(5.38) 



dz 



which shows that the boundary condition (5.6) is satisfied. Analogous considerations can be done for 
proving (5.7). Theorem |5.5] is then proved. □ 

Proof of Theorem 5.6. Following the steps in ]I(J we consider the polynomial 

Z g {\) := -xP 9 (X) - \2tP 9 (A) + R g (X), 
where R°(X) — f(u) and R 9 {X), g > 0, is given by the expression 



(5.39) 



2g+l 2g+l g k 

R^X) = 2j29u k q g (n) ]J (A - u,) +^? fe (fi)^(2Z - l)f fc _,i* X (A), (5.40) 

k=l l=l,ljtk fc=l Z=l 



with the polynomials P ; S (A), Z > 0, defined in (2.8) and the functions qk{u), k = 1,. .. ,g, defined in 
(5.8). Then the hodograph transformation (5.3) is equivalent, for g > 0, to the equation 

Z 9 (X) = 0, g > 0. (5.41) 



The proof of the above proposition is obtained observing that the u>i(it)'s defined in (5.E) are given by 

RHui) 

the ratio Wi(u) — -^r, — r, % = It 



,2(7 + 1, where R 9 (X) is the polynomial defined in (5.40). Hence 



we can write the hodograph transformation (5.3) in the form 

[-xP 9 (X) - UtP? (A) + R 9 (X)] x=m =0, i = 1, . . . , 2g + 1. 



(5.42) 



For 5 > 0, Z 9 (X) is a polynomial of degree 2<? and because of ( 5.42 ) it must have at least 2g + 1 real 
zeros. Therefore it is identically zero. Putting equal to zero the first g+1 coefficients of the polynomial 
Z 9 {X) and using repeatedly the identity (5.13) we obtain ( 5.15| ). 

Putting to zero the coefficients of Z 9 (X) of degree to degree g — 1 is equivalent to the equations 



Z 9 (X) 

ak w(A) 



dX = 0, fc = 1, 



(5.43) 



23 



Indeed when (15.15) is satisfied, the differential — — — dX is an holomorphic differential and by (5.43) 

y(x) 



it has all its alpha-periods equal to zero. Therefore it is identically zero. Hence the last g equations 
reads 



-xP${X) - 12iPf (A) + RB{X) 



a k 2/(A) 
= / 2 



y(A) 

ezv 3,%( "in^/i(A-«i) 



(5.44) 



dA, fe = l, ...,g. 



In the third equality of ( 5.44 ) we have used the fact that 

if (A). 



/ ^vdX= [ <jf(A)=0, l>0, fc = l,...,ff, 



because of the normalization conditions (2.6). The function 5 ,ff (A;u) defined in (5.19) satisfies the 
relations 



A 



A - 



= 2G) Ui * ff (A; 3), 2d Ui q g (u) = u). 



(5.45) 



Using (5.45) we can rewrite the last term in (5.44) in the form 



0= 2- 



dX, k = l,...,g, 



=2 



«2fc-l 



«2k-l 



4/ y(A) 9 A *f(A;3) + ^9 Ui *ff(A;u) ) dX, k = l,...,g, 



"2k 



where the last equality has been obtained integrating by parts. Using the definition of $ 9 (A;u) in 



(5.17) we rewrite the above relation in the form 



= -4/ y(X)^(X;u)dX, k = l,...,g, 

'«2Js 



which is equivalent to (5.16) 



(5.46) 

□ 



5.1 The function Q' and linear-overdetermined systems of Euler Poisson 
Darboux type. 



In this section we show that the set of algebraic equations described by the moment conditions (4.17) 



and the normalization conditions ( 4.13 ) can be written in terms of solutions of linear overdetermincd 
system of Euler-Poisson-Darboux type introduced in pT| . 
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Theorem 5.7 For any monotonically increasing analytic initial data satisfying (3.1), the set of alge- 
braic equations defined by the moment conditions (l-l r i) is equivalent, for g > 0, to (5.11); the set of 
algebraic equations defined by the normalization conditions ( J^.lh ) is equivalent, for g > 0, to ( 5.1i ), 
namely 

2g+l 



Y] 9 v>j<l9-k{ti) - kq g - k +i(u) = — I 



(p'(z) - a\z))z k 



2g+l 



Y] d u qi(u) - (g - l)q 2 (u) - 6< = — 

3=1 

2ff+l 2 S +1 

2 Ujd Uj qi(u) + q^u) — x — 6t y~] Uj = — 

3=1 m 



(p'(z)-a'(z))z^ 



dz, k = 0, . . . , g — 2 



!/+(*) 



(5.47) 



«2fc-l 



(p'(z)-g'(z))z* 



2y + (A)$(A;u)dA- / (G' + (X) — Q'_(X))dX, k = 1, 



Proof: we first consider the moment conditions. Since a'(z) = 6tz? + — z 2 we have 



a'{z)z l 

y( z ) 



dz = 0, k = 0,...g-2 



1 

7TZ 



— rfz = 6i, — 

y^yz) tti 



a'(z)z & 



2g+l 



(5.48) 



dz = 



6*E 



For analytic initial data we can deform the integrals 



p'(z)z k 



p'(z)z h 



dz, fc = 0, . . . , g 



(5.49) 



where So is a close loop around the interval [0, "29+1) passing through zero and the dj are close loops 
around the intervals Ij and within the domain of analyticity of the initial data. We can deform the 
contours Udj to a single contour C = U^ =0 o7j U| =1 jj where the contours jj are the closed contours 
plotted in the figure below. 




a„ 



The integrals 



therefore 



p'{z)z h 



dz = 0, j = 1, 



p'(z)z k 1 f p'(z)z k 



. -dz 



-dz, k = 0, 



(5.50) 
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dz, 



I -n p' ' (z)Z 

From the right hand side of ( 5.5C ) it is straightforward to verify that the integrals J c — 

A; > 0, are symmetric with respect to the variables u\, . . . , W2n+i and satisfy the linear overdetermincd 
system of Euler-Poisson-Darboux type defined in (5.8) with the initial data 

P ' {z)zk i Z 

[«l=«2 = ---=«2g+l=«] 

= 2ud u F 1 (u) + F 1 (u), 



1 



2 

m 



p'(z)z 9 



d u F q _ k {u) - fcF„_ fe+ i(u) 



dz 



(5.51) 



[Ul=U2 = 



="2 g + l="] 



where the functions Fk(u), k = l,...,g, have been defined in (5.8). Because of the uniqueness of the 
solution of the boundary value problem (5.8) the first g + 1 equations in ( |5.47 ) are satisfied. 

As regarding the normalization conditions (4.13) we have the following identity for A € Ij, j = 
0,...,g 



g' + (x)-g'_(x) = 



2?/ + (A) 



v. p. 



y + W 



P'(z) - a'(z) 
(z-X)y+(z) 

^Jd } {i-X)y+(z) 



dz 



dz 



(5.52) 



y+(X) fp'(z)-a'(z) 



ni J c (z - X)y+(z) 
2y + (\)(<S> g (\;u)-6te 



dz 



g0), 



where the function Q 9 (X;u) has been defined in ( 5.17 ) and e g o is equal to one for g = and zero 
otherwise. The last identity in ( |5.52| ) has been obtained performing computations similar to the ones 
in (5.51). Therefore 



0= / (S;(A)-Sl(A)dA = 2 



y + (X)&(X;u)dX, j 1 «/. g -0 



and theorem 5/7 in then proved. 

Theorem 5.8 The relations ( 5.J^\ ) are satisfied for smooth initial data. 



Proof: the proof is obtained combining theorems 4.3, 5.5, 5.6 and 5.7. 



(5.53) 

□ 

□ 



In the following we write the variational conditions in (3.7) and (B.8) in terms of the function $ 9 (A; u) 
defined in (5.17). 



Theorem 5.9 The variational conditions (3.7) and (S.t) can be written in the form 

g + (A)-a_(A) 



> 



2i 



zT(£)($ 9 (£;w)-6ie g0 R, A e (u 2j , u 2 j-i), j = l,...,g + l (5.54) 
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and 



0<g' + +g'_-2p' + 2a' = 2y(A)($ s (A;u)-6te s0 ), Ae {u 2j+l , u 2j ), j = 0, . . . ,g, u = 1, (5.55) 



where the function $ 9 (A;u) /ias 6een defined in ( 5.11 ) and e g o is equal to 1 for g = and zero 



Proof: we first prove (5.54). From (5.52) and lemma 4.1 we obtain 



G+-G- 
2i 



g> + -g>_ 
2i 



,j y(0* 3 foiZ)<%<0, A G yj 9 j t\{u 2 j 1 u 2j -i) 



which coincides with (5.54). 
As regarding (5.55) we have 

2p'(z) 



G' + + GL = y -^[ 
m Jc, 



y+(z)(z-X) 



dz - 2a' (A) - 12te g0 y(\), \€ (u 2 j+x,u 2 j), j = 0,...,g, (5.56) 



where e ff o is equal to one for g = and zero otherwise. 
In order to express the integral 

2p'(z) 



-dz 



y+(z)(z-X) 

as the solution of a linear overdetermined system of Euler-Poisson-Darboux type, let us consider the 
function g g+ i(ui, u 2 , u±, u 2 , . . . , u 2g+ i) = g s +i(iti, u 2 , u) which satisfies (5.8) with initial data 
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q g+1 (u,u, ...,u) 



2g+3 



(25 + l)H 



f {9) (u). 



Then 



'2g+l 



^ 3 Mj q s+ i(ui, M~2, u) + ^Og+^ui, u 2 , u) + du 2 q g+1 (u 1: u 2 , u) 

3 = 1 



<f> 9 (\;u). (5.57) 



The derivation of the above identity is straightforward. Next we consider the first identity in (5.47) 
for the function g g +i(ui, u 2l u). Let us suppose u 2 j > u\ > u 2 > u 2 j+i, for j = 0, . . . ,g, uq = 1, let us 
define the interval / = (u 2 ,u~i) and Y(X) = y{\)^ (z — u 2 )(z — u~i). The function Y(\) is analytic in 
the complement of (— oo, 0] U C g U / and positive for A > u\. Y + (X) denotes the boundary value from 
above. Then by ( [5.47 ) we have 



2g+i 1 
2] d Uj q g+ i{ui 1 u 2 ,u) + d^qg+iiui, u 2 , u) + du 2 q g+ i(ui,u 2 , u) = — I 
j= i m J c a iii 



Y+(z) 



dz. 
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The limit Ux — > u 2 of the left hand side of the above identity has been obtained in (5.57). Here we 
derive the limit of the right hand side. Let us define u 1 = A + VS, 1I2 = A — \f&. Then 



lim — 

(5^0 7TI 



c a ui 



Y+(z) 



1 



w d.+iimi r vo ^i dz 



2p'(z) 



Combining the previous three relations we obtain 



dz — 



V(A) 



, A e (ti2i+i,u 2j ), j = 0, . ..,5, u = 1. 



(5.58) 

□ 



£s y+(z)(z-A) y(X) 

Combining ( 5.56| ) and fl5.5Sj ) the variational condition (5.55) can be easily obtained 

Remark 5.10 We observe that the function $ 9 (A;u) — 6tS g o must be nonzero at the branch points. 
This is obvious for g=0. Indeed $°(u;u) — 6t = f'(u) — 6t = l/u x (x,t). For g > the derivatives 
d x Ui(x,t) are given by the expression J7fi|/ 



d x Ui(x,t) 



2Y\fJ[\u l -u 3 )^^ Ul -u)' 



i= 1,...,2$ + 1, g>0, 



where Pq has been defined in (2.8). From the above expression it is obvious that the d x Ui(x,t) 's are 
nonsingular if <& 9 (ui] u) ^ 0, i = 1, . . . , 2g + 1. 



5.1.1 Phase transitions 

We derive the equations which determine a change of genus of the solution of the hodograph transfor- 



mation (5.3) or the set of algebraic equations (4.17) and (4.13). We observe that a transition necessarily 



occurs when one of the two conditions ( 5.54 ) or ( 5.55 ) fail to be satisfied. If the conditions ( 5.54 ) fail 
to be satisfied at some point v in the bands (u2j, U2j-i), j = 1, • • • , g + 1, it follows that 

23+1 y(0(® g (Z;u)-6te g o)dt; = o, 

(5.59) 

$ 9 (v;u) - 6te g0 = v € (u 2 j,u 2 j-i), j = 1, . . . ,g + 1. 
Indeed v must be a stationary point of the integral J^ 2j+1 ^0 — 6ie 9 o)rf£- When we can solve 



the above system together with the moment conditions ( 4.17 ) and the normalization conditions ( 4.13 ) 
for some t > 0, we obtain a point x(t), v(t) and u\{t) > u 2 (t) > • ■ • > U2 g +\{t) of the boundary 
between the g-phase solution and the (g + 1)— phase solution. 



In the same way if the conditions (5.55) fail to be satisfied at some point v in the gaps (v,2j, u 2 j+i), j 
0, . . . , g, it follows that 

$ 9 (u;u) - 6te g0 = 
d v <f> 9 (v;u) = 0. 



(5.60) 



The above system describes the point of phase transition when a new band is opening. 



Systems (5.59) and (5.60) have been obtained in |10[ studying directly the hodograph transformation 



5.2 Lax-Levermore-Venakides functional 

Following the steps in |jl2"f we construct the maximizer of the Lax-Levermore-Venakides functional 
Let us make the change of variable A = 1 — rj 2 , taking the upper complex half plane onto 
C\(— oo, 1]. The function Q(X) transforms to ^(rj) = <?(A). The function T is analytic off the real r\ 
axis. We extend the definition of T onto the lower complex 7/-plane by the relation T{—r}) = —T{r\). 
The RH problem for the function T is T+ — T- = G+ + G- and T+ + T- = — (G+ — G-)signi]. On 
the real axis we define the function tp*(n) to be equal to zero outside the interval [—1,1], while on the 
interval (—1, 1) 

1 



The function ij)*{rf) is the unique maximizer of the functional B], [jl2l 



(5.61) 



Q(V0 = -[(2a, V) + (14, i>)], V € ^([0, 1]), V < 0, 

7T 



where 



i 

n 





7?-/i 


/o 





-0(m)o> 



0(77, x, i) = 4t?7 3 — xr\ — Qtrj + 



1 



/(0 



1— 77- 



1 



The maximization problem is attacked analytically by solving the variational conditions 

Lip(fj) + a{f], x, t) = 0, when ipif], x,t) < 
Lip(rj) + a(r), x, t) > 0, when ip(ri, x, t) = 0. 



(5.62) 



Namely if tp satisfies (5.62), then ip = tjj* [J4|. In [|12| it is shown that the variational conditions (5.62) 
transform to (3/7) and (^8). The benefit of the variational formulation is that, due to the convexity 



of the maximization problem, the uniqueness of ip* and hence of G is guaranteed. Furthermore for 
each fixed x and t the support of the maximizer in uniquely defined. An important consequence of 
this result is the following. 

Theorem 5.11 For almost all x and t > the solution u\(x, t) > U2(x, t) > ■ ■ ■ > U2 g +i(x, t) of the 



hodograph transformation (5.5) 

x = —Vi(u) + Wi(u), i = 1, . . . , 2g + 1, g > 0, 

which satisfies the constraints (3.7) and (S.h ) or, equivalently, ( 5.54 ) an d A?.<?4 A exists for some g > 
and it is unique. 



29 



As a final result we write the Lax-Levermore-Venakides maximizer as the solution of a linear overde- 
termined system of Euler-Poisson-Darboux type. 



Theorem 5.12 For g > the maximizer ip*(rf) in (5.61) can be written in the form 

flfc+i) = 0, ne (0i,/3 2 )U(/3 3 ,/?3)U---U(/3 2ff+ i,l) 
r(m ft, ft, ■ ■ .,p2g+i) = -2r,&(l - V 2 ; 1 - 0{, 1 - . . . , 1 - /3| 3+1 )y(l - V 2 ) 
when i] £ [0,/3i] Uj =1 [/%, /%+i], 

/3 fe = VI -u k , k = l,...,2g + l, 

2g+l 
3 = 1 



and the function $ 9 has been defined in ( 5.17 ) 



We observe that the condition ( 5.55| ) implies ip*(r]; 01,02, ftg+i) < for rye (0, /3i)U| =1 (/3 2 j, Ay+i) 



and vice-versa. The proof of the theorem follows directly from (5.55) and ( |5.6l| ) 



Remark 5.13 All the information on the initial data of the Lax-Levermore maximizer is contained 
in the the function <5> ff (A;u). 

In |7Z|/,.J7^/ an upper bound to the genus of the solution of the Whitham equations was provided. We 
selected initial data such that the maximum number of real zeros of the function $ s (A;u) — 6tS g o is 
2N — g for < g < N where N is the supposed upper bound of the genus. Then we showed that, in 
such a situation, phase transitions from solutions of genus g < N to solutions of genus g > N do not 
occur. The proof of the theorem is obtained only studying the Whitham equations and the hodograph 



transformation (5.5 ). 

In this new representation of the Lax-Levermore- Venakides maximizer, the above result has an easy in- 
terpretation. Namely it gives an upper bound to number of intervals where the function y(\)((fr 9 (A; u) — 
6t<5 s o) is positive for < g < N . 



6 Conclusion 

In this paper we have studied the Cauchy problem for the KdV equation with small dispersion and 
with monotonically increasing initial data. We have used the formulation of the Cauchy problem for 
KdV as a Riemann-Hilbert (RH) problem Jl8| . The small e-asymptotics is obtained using the steepest 
descent method for oscillatory Riemann-Hilbert problems introduced in [jl6| . Analyticity of the initial 
data is essential for this step. 

This procedure leads to a scalar RH problem for a function Q and a set of algebraic equations con- 
strained by algebraic inequalities. 
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The scalar RH problem for the function Q is well defined even for smooth initial data. 
In this paper we have shown that, for smooth monotonically increasing initial data bounded at infin- 
ity, the set of algebraic equations obtained through the Deift, Venakides and Zhou approach can be 
expressed as the solution of a linear overdetermined systems of equations of Euler-Poisson-Darboux 
type. We have also shown that this set of algebraic equation is equivalent to the set of algebraic 



equations defined by the hodograph transformation (5.3). The uniqueness of the Lax-Levermore max- 
imization problems guarantees the uniqueness of the solution of the Cauchy problem for the Whitham 
equations. As a final result we have shown that the Lax-Levermore maximizer can also be expressed 
as a solution of a linear overdetermined system of Euler-Poisson-Darboux type. 
We believe that the above results can be extended to bump-like initial data. 
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for oscillatory Riemann-Hilbert problems. 
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